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THE MOTION OF A SATELLITE. 



BY PEOF. ASAPH HALL, NAVAL OBSEEVATOEY, WASHINGTON, T>. C. 

(1). The apparent orbit of a satellite around its primary will be in 
general an ellipse, which will be the projection of the real orbit of the sat- 
ellite on a plane perpendicular to the line of sight, or to a right line drawn 
from the observer to the centre of the planet. Assuming that the orbit of 
the satellite around its primary is known, and that its motion is purely el- 
liptical, or free of perturbation, we have the problem of determining from 
the known orbital motion of the satellite its apparent motion. The solution 
of this problem has been given by Bessel in the Astronomische Nachrich- 
ten, Bd. 9, p. 7; and again in his Astronomische TJntersuchungen, Bd. 2, p. 
37; and also by Marth in the Astronomische Nachrichten, Bd. 44, p. 113. 

Denote by s and p the apparent distance and angle of position of a sat- 
ellite, or its polar coordinates referred to the centre of the planet. The 
zero of the angle p is the circle of declination passing through the centre of 
the planet, this angle being counted from the north toward the east, and 
from 0° to 360°. 

Let a and d be the right ascension and declination of the planet, a' and 8' 
the same quantities for the satellite. From the spherical triangle between 
the pole of the equator, the planet and the satellite, we have 

coss = sin 8sin ^'-j-cos ^cos 8' cos (a' — a) ^ 

sin s cos p = cos 8 sin 8' — sin 5 cos 8' cos (a ' — «) [■ . (1) 

sinssinp= cos 5' sin (a' — a) J 

If we denote by r the radius vector of the satellite, and by a and d its 
right ascension and declination seen from the planet, and also by p and p' 
the distances of the planet and the satellite from the Earth, we shall have 
by equating the values of the rectangular coordinates. 
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p' cos B' cos a' = pcos§ cos a -\- rcosd cos a, 

p' cos d' sin a' =■ p cos 5 sin a -|- r cos d sin a, 

p' sin ^' = /) sin 3 -f »• sin d, 

or by a simple transformation 

p' cos d' cos («'— a) = /) cos 5 -|- r cos d cos (a — a), 

jo' cos 5' sin {a! — «) = r cos d sin (o— a), 

p' sin ^' = /? sin ^ + r sin d. 

If we multiply equations (1) by ^' we have by means of the last equa- 
tions 

/>' cos s = |0-fr . [sin d sin ^-f-cos d cos d cos(a — a)] "^ 

jo' sins cos p= r. [sin d cos 5 — cos d sin 5 cos(a — a)] )-. (2) 
p' sin s sinp = r . cos d sin (a — a) J 

If we place now 

f = - . cos d sin (a — a), 

3y = - . [sin d cos <? — cos d sin ^ cos [a — a)], 



(^ = - . [sin d sin 5-fcos d cos 5 cos (a — a)], 



equations (2) become 



We have therefore 



p'&^ss = p.(l-\-C) ] 
^' sins cos p = />.;y )■ . (3) 



p' sin s sin p = />. ? J 



tang s sin » = — - — , 
tang s cos p = 



1 + C 

If we assume tang s = s, which may generally be done without sensible 
error, and put a; = ssinp, y = s cos p; we shall have in seconds of arc 

1+C ^ 1+C 

where to .-= 206264.8. The equations (3) give the values of p', s and p 
when those off, tj and C are known. It is evident also that C has no influ- 
ence on the value of p, and since it is generally very small with respect to 
unity its influence on the value of s is slight. 
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We have now to express the quantities f, -fj, <r, in terms of the elements 
of the orbit of the satellite. Continuing the use of Bessel's notation, let N 
be the longitude of the ascending node of the orbit of the satellite on the 
Equator, and / its inclination to the Equator. Let u be the argument of 
the latitude, or the angle between the node and the satellite counted in the 
direction of the motion. From the right angled spherical ti'iangle having 
M for a hypotenuse, and the sides d, and a — iV, we have 
cos d sin (a — N^ = sin u cos J, 
cos d cos {a — N) = cos w, 
sin d = sin u sin J. 

Eliminating N from the left hand side of these equations we have 
cos d cos {a — «) = sin u sin (« — iV) cos /-fcos u cos {« — N), 
cos d sin (a — «) = sin m cos (« — iV")cos / — cos u sin (a — iV), 
sin d = sin M sin J. 

Denote the mean distance of the planet from the sun by jOq, and the semi- 
major axis of ths orbit of the satellite by ^Oj sin //; if we express the radius 
vector r by means of the semi-major axis the value of r will be 

jOq sin A .r, 

or approximately p^ 



.r. 



The values of f, t], ^, are therefore 

w.f = r .f^A. -^ sin mcos (« — N ) cos J — cos u sin (« — N) > 

w.fj = r . C-O- J. ^ sin M [cos d sin J — sin d cos Jsin (« — N )] 



}. 



(4) 



— cos u sin d cos (a — N) )■ 

oi.!^ = r.^A. -{smu [sin d sin /-f-cos d cos /sin (« — iV)] 

f^ -\- cos U COS <? COS (« — N) \ _ 

Introducing the auxiliary quantities /, JF, g, Q, h, H, by the formulae 
sin /cos F = cos (« — N) cos J, 
sin /sin i'' = — sin (a — iV), 
cos / = — cos (a — N) sin J", 

sin g cos G* = cos d sin /—sin <5 cos /sin (« — iV), 
sin gsmG — — sin d cos (« — iV"), 
cos g = cos 5 cos /-|-sin ^ sin /sin (a — iV"), 

sin hcosH= sin 5 sin /+cos d cos /sin (a — N ), 
sin A sin jy = cos ^ cos (a — iV^), 
cos h = sin 5 cos / — cos 8 sin /sin (a — iV), 
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we have 



a).i = r. £fi Jsin/sin (F + u) 
(ti.Tj = r. lO dsmg sin{0 + u) 



V. 



(5) 



w.r = r .^ A sm hsiQ. (H 4- u) 
P 

Denoting the mean longitude at the instant from which the time t is 
counted by E, the daily motion by k, the longitude of the perihelion by P, 
the eccentricity of the orbit by e, the mean, the eccentric and the true anom- 
alies by fi, e and f, we have 

E-\-M — P = e — esins ^ 



f^ = 
tang 



r = 1 — 



+e 



—e 

e cos £. 



. tang j£ y . 



(6) 



The preceding formulae are those given by Bessel. The six auxiliary 
quantities /, F, g, Q, h, H, are similar to those introduced by Gauss for 
computing the place of a planet (Zach's Monatliche Correspondenz, Mai, 
1804, or Werke, Bd. 6, p. 94), and their geometrical meaning is readily 
s een. Thus, if we take the origin of rectangular coordinates at the centre of 
the planet, the line drawn from the planet to the Earth as the axis of x and 
take the axis of y in the plane of the declination circle, we shall find that 
/. g, h, are the arcs drawn from the north pole of the orbit of the satellite 
to the points where, the axes of x, y and z cut the sphere; and that F, O, 
H, are the angles at this pole between the arcs drawn to the north pole of 
the Equator, and those drawn to the points where the sphere is cut by the 
axes of X, y and s. 

By introducing a single additional auxiliary quantity Marth has given a 
very simple and elegant method of computing s and p. Thus if we com- 
pute the angle h by the formulae 

sin A sin ^ = cos (a — N) sin J 
sin A cos ^ = — sin (a — N) sin J sin d — cos J^cos d . 
and denote by a the supplement of the angle between the line drawn from 
the Earth to the planet and the radius vector of the satellite, we have from 
the spherical triangle formed by the directions of the selines and the pole 
of the satellite's orbit, 

sin <r sin {p — h) = cos (H-\-u) s. 

sin a cos (p — h) = sin {H-\-u) cosh V . (8) 

cos <T = sin (if+w) sin h * 



}> 



(7) 



—133— 



We evidently have 



p' sin s = r sin a, 
p' cos s = r cos <T-{-p, 
and therefore we can find p', s and p from the equations 

p' sin s sin (p — h) = r cos {II-\-u) •\ 

jo' sins cos (p — i) = r sin (iJ-fw) cos A V. (9) 

p' coss = r sin {H-\-u) sin h-{- p ^ 

By putting 

sin (jET+m) sin A-|-2 = - 
r T 

and eliminating p' we have for computing s and p, 
tang s sin (p — h) = r cos {II-{-u) 
tang s cos (p — k) = r sin (H+m) cos A. 
If we wish to use the angle ff we may compute s by means of the equa- 

tions sin s = - .sin (<t — s), 



and for p' we have 



P 



sni (T 



sin (<T — s) 
Both these equations come directly from the plane triangle between the 
Earth, the centre of the planet and the position of the satellite in its orbit. 
The preceding formulae given by Marth are rigorous, but generally it will 
be sufficient to use the first two equations of (9). If we neglect the diifer- 
ence between p and p', and express s and »• in seconds of arc, we have 



s sin (p — h) = - . cos {H-{-u) 
P 

s cos (p — Ic) = ~ . sin (If+tt) cos h 
P 



(10) 



Equations (10) are convenient for use if from the observed values of s 
and p we wish to compute r and u. For this purpose we may use also the 
auxiliary quantities of Bessel, which give 

r . sin il{F+ (?)+«] = % . sing; sin p+sin/cos_p 
L2\ -r ;-r J 2 sin/sin^rcosKJ'— (?)' 

r . cos [i(-F+ G-)+u-] - ' . sinffsinp-sin/cos^ 
L2V -r ;-r J ^ sin/sin ^ sin ^(i''— (?) ' 

but it is more convenient to use the formula of Marth. 

In computing the quantities /, F, g, G, h, H, we use the quantities J 
and N, which are themselves subject to slow changes on account of preces- 
sion. But we may easily form a table of the values of J" and iV which will 



—134— 

enable us to take account of these changes in the values of the auxiliaiy 
quantities. If moreover we wish to take account also of the nutation, and 
denote by Jco and Jl the nutation in the obliquity of the ecliptic and the 
longitude, we may employ the differential equations of a spherical triangle. 

[f u) be the obliquity of the ecliptic we have 

dJ = cos N. A(D — sin iVsin w. Jl, 
J sin N . Jo) 4- cos iVsin co . Jl 

sin J 
dN = — cos J. dw -{-cos a) . Jl, 
or dN = — (sin -^7'. Jco-}-cos Nsin co . Jl) cotg /+cos co .JL 

The first and last of these formulae will give the corrections to the values 
of /andiV. 

(2). If we assume values of the elements of the orbit of a satellite, the 
preceding formulse enable us to compare our observations with the elements. 
In order to correct the assumed elements we must form equations of condi- 
tion by differentiating the first two of equations (5). Denoting the partial 
differential coefficients by a, b, c, d, e, f; and by v the residual = compu- 
ted — observed value, the equation of condition will be of the form 
a . dE+b . edP-f c . de-\-d.dJ-\-e . sin JdN+f.dJ+v — 0. 

Each observation gives two equations, one from the value of ?, and one 
from the value of t^. The changes in the value of ^ can be neglected. 

If we denote by n the longitude of the ascending node of the orbit of the 
satellite on the Ecliptic, by i its inclination to the Ecliptic, by w the 
distance of the Ecliptic from the Equator, measured on the orbit of the 
satellite, by v the longitude of the satellite in its orbit, we have 

V — n-f-w = u, and f = v — P. 
If we put 

f = PAJsmf; q' = P^Jsmg; h' = P±Jsmh, 
p P P 

F'==F-\-yv—n; G'=- G-{-w—n; H'= H-\-w—n, 
the values of f, yj and ^ become 

co.$ = r.f' sin {F'-^-v), 
CO .7j ^= r .g' sm{G'-j-v), 
(0 . !^ = r . h' sin {H' +v). 
The value of « . f may be written in the form 

CO . ? -/'. [sin {F'+F). r cos (v— P)+cos {F' + P) . r sin («— P)]. 
Differentiating with respect to the elements in the plane of the orbit we have 
w.d^=f sin {F'+P) . d[rcos(v— P)]+/' cos (P'+P) . d[r sm{v—P)'] 

+rf' cos {F'-l-v).dP. 
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We have also, since v — P is the true anomaly, 

r . cos {v — P) = cos £ — e, 

r . sin {v — P) = sin e V(l — e^}, 

E — P = £ — e sin £, 
and therefore 

d . [r cos (» — P)] = — sin £ . ds — de, 

d . \r sin (»— P)] = cos £ . ^(1— «') • dz-^±^, 

d . (E- — P) = (1 — e cos e) . & — sin £ . de, 
or 

de = d{E—P)-}-Buiz.de 
1 — e cos £ 
Eliminating ds we have 

d . [^ cos (^;-P)] = -— ?^"^^— . d (jB;— P)— (l +_.^yif!— ). de, 
1 — ecose V 1 — ecos£/ ' 

d . [r sin {v-P)-] = coss- V(l-e^) ,a{E-P) 
1 — ecose 

. /cos£. |/(1 — e^) e \ . , 

+ — 1 — — ^— 77T — 2\ • sin e de. 

\ 1 — ecos£ |/(1 — e^)/ 

But since 

_£."_i_ = ^^^#=::?^ : and cos {v-P) = ^^^', 
1— ecos£ i/(l — e') 1 — ecos£ 

the values of these differentials become 

d . Ir cos (._P)]= -^i^^) . d iE-P)-[l+^^^^:zp).de, 

d.[r sin iv-P)-]= gg^^H:-^ ■ d(^-P)+ "°;;(';1^--^^ .de. 

Multiplying the first differential by /' sin (P'+P), and the second by 
/' cos (P'4-P), the value of w .d$ becomes 

(o.di = ^^^ . I cos {F'+v)+e cos (P'+P) I . dE 

+ ^7?iz:i2-) • { <^0'^ (^'+^) ['■^'(l -e2)-l]-e cos (P'+ P) | . rfP 

+ 77r=^2y {cos(P'+t)) sin £ — sin {F'+P).^{l—e')\.de. 

If we substitute the value of r = 1— ecos £, in the coefficient of dP, and 
observe that the value of the partial differential coefficient d is evident, we 
shall have the values of the first four coefficients given by Bessel ; Astrono- 
mische Nachrichten, Bd. 9, p. 10. 
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In order to find the two other coefficients, we take the vahie of a> . f 

CO .^ = r .^ J .[sin u cos (a — N) cos J — cos m sin (a — Nj]. 

Differentiating with respect to the elements that determine the position 
of the plane of the orbit we have 

o) .d$ = r.^J.< [cos «cos (a — N) cos /-|-sin u sin (a — NJ] . du 

-\- [sin u sin(a-iV)cos J-|-cos u cos(a-iV)]diV-sin u cos{a-N) sinJdJ J- 
The values of the auxiliary quantities / and i^give 

o) .d? = r.^A . i sin/cos {F'-+v). du+sinfcos {F' +v) cos J . dN 
P ^ 

— cos u COS /sin J . dN-\-sm u cos/ . dJ > . 

If we now consider the spherical triangle formed by the equator, the 
ecliptic and the orbit of the satellite, we shall have by the formulae for the 
differentials of a spherical triangle, 

dw = cos i . dti — cos J . dN, 

dn = -'^. dJ +S^Z^ . dN, 
sin I sm I 

J sin w cos i , 7^ . cos w cos t sin / , ^r r j at 

or dw = -. — -. — ■. dJ + i — . ^dN — cos J. dN. 

sin z sin t 

We have therefore 

d (w — n) ^ ( — cos J — cos w sin J tang Ji) . diV^-t-sin w tang ^i . dJ. 
Since du^= d (w — n), by substituting this value of du we have 

o) . dq— — sin/. < /'.tang Jicosw.rcos(J^'+w)+L_o Jcos/.rcosit > . dN 

+ \ /'. tang Ji sin w,*" cos (i^'+'y)H-^z/ cos /.r sin tt \.dJ. 

The values therefore of the partial differential coefficients a, h, c, d, e,f, 
are as follows : 

a = 0.. 4t = TTT^^-I «o« {F'+v)-\-ecoB{F' + P) | 



dE ViX—^^) 

_ d^ _ -f 



b = M.—^r^ 



edP i/(l — e') 



,| cos(i^'+«;Ucosev/(l— e^) 



l+v'(l-e^) 

. ^ r>ns ( VA-ii 

de V{l—e') 



TwJ+-(^'+^) 



:a>.-^= ,{' .„ . / COS (F'+v) sine— sin (i?''+P) • t/(1— e' 
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da A ' 

e = CO,— — TyTT= — /'. tang|icosw.rcos(i'''+'i;) — C-& J. cos/. rcosM 
SID odJy p 

f = 0). -jj = /'. tang Ji sin w.r cos (i'''+»)+—^- cos/, r sin M. 

(11) 

In order to obtain the coefficients in 5f we have only to put g' and G' in 
the place of /' and F\ 

The preceding formulEe give the values of the differential coefficients 
adapted to the quantites 

fr=a; = ssia^: y^ = y = scosp; 

and we must therefore combine our observed quantities s and p to find the 
residuals for the equations of condition. Generally it is better to compare 
directly with each of the observed quantities, since one of them may be er- 
roneous and the other not. If we wish to do this in the present case, we 
must use the formulae 

ds = sinp.dx + cosp . dy, 
s, dp = cos p. da; — sinp ..dy. 
By means of equations (1 1) the value of dx is 

dx = a .dU + b .edP + o .de + d.dJ + e.sinJdN + f. dJ, 
with a similar value for dy. Substituting these values for dx and dy in the 
preceding equations we shall have the equations of condition for ds and 
sdp. If we make 

m sin M = f' cos {F'-\-v), n'sm N' =/'cos {F'-\-P), 
m cos M = g' cos (G' -\-v), n'cos N' = g^'cos {G'-{-P); 

I sin L =/'sin (F'+P), k'sin K' = P^A cos/ 

IcosL = g' sin ( G'+P), ¥ cos K' = P^A cosg. 

The values of the differential coeflScients for ds and sdp are as follows, 
where sin f' is put for e on the right-hand side of the equations. 
For the distances ; 

di ^ -^^ 7.|mcos(lf— p)+w'cos(iV'~p)sin^'| 

-r^ = — ■< m cos {M—p) (cos e cos y'+tang ^f') +n'cos {N' — p) > 

~r- = y A m COS (M — p) sin e — I cos iL — «) cos w' \ 

de cos f' { \ »-/ \ JTJ T j 
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ds s 

ds 



sin JdN 
ds 






dJ ^' 



'. < m cos (M—p) cos w tang |i -j-k' cos (^' — p) cos m >■ 
■J m cos {M—p) sin w tang |i+A' cos (^' — p) sin w } • 



For the angles : 

sdp 

dE cos 



(12) 



; — J A m sin (if— p)+n' sin {N'—p) sin f ' |- 
iSP "" ^^V • { '^ s^° (lf~p)(cos e cos^'+tangi9j')+»i'sin(iV''— p) | 

-J- = ^ . K m sin (if— p) sin e — ^sin {L — p) cos^' J- 

sdp 



dJ 
sdp 







■ Tdw'^ — ^ ') "'' ®^° {^~~P) cos w tang |i + Z;' sin (JST' — p) cos « > 

-_-E = r • ^ m sin (if — p) sin w tang Ji + k' sin ( JT' — p) sin M >■ 

(13) 
Generally in equations (12) and (13) we may put <p'=0, and these equa- 
tions will become simpler. 

The values of i, n and w will be found from the spherical triangle 
between the Equator, the Ecliptic and the orbit of the satellite. This tri- 
angle gives by means of the Gaussian equations, if we denote by o) the ob- 
liquity of the ecliptic ; 

cos |i sin |(w-f n) = — cos JiVcos |(<u — J) ^ 
cos J« cos J(w-f w) = sin ^^ cos |(tu+J) ' 
sin |i sin J(w — n) = cos JiV^sin |(a> — J\ i 
sin |icos |(w — n) = — sin JiVsin f(o>-f J) J 

(3) In what precedes it is assumed that the elements of the orbit of the 
satellite are known, and the formulae given are sufficient for comparing the 
observed positions of the satellite with those computed from the elements, 
and for forming the equations of condition necessary to correct the assumed 
values of these elements. In order to obtain an approximate knowledge of 
the orbit of a satellite various methods may be used, but in this case, the 
elaborate formulae given for computing the orbit of a planet are not neces- 
sary ; and in my own experience I have found it nearly useless to employ 
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them. The reason of this is that the small errors of observation exert too 
great an influence on the resulting elements, so that it is better to use sim- 
pler methods. If we assume the orbit of the satellite to be circular, and 
d(itermine by observation the form of the apparent ellipse described by the 
satellite, and the position of its major axis, the solution of a spherical trian- 
gle will give the values of N and J, the quantities which fix the position of 
the plane of the orbit. With these quantities we can compute the values 
of the auxiliaries /, g, h, &c., and by means of equations (10), and the 
observed values of s and p, we shall have the values of r and u, the polar 
coordinates of the satellite in the plane of its orbit. It is then easy to de- 
termine the elements of the ellipse described by the satellite around its pri- 
mary, the position of the line of apsides, and the place of the satellite in its 
orbit at a given time. 

The following figure explains the geometrical meaning of the auxiliary 
quantities ; and one will be able to deduce the formula for computing them 
directly from the spherical triargles. 

TZ=f, 

TY=g, 

TX=h, 
OTZ= F, 
OTY= G, 
OTX=H; 

0T= /= in- 
clination of orbit 
to Equator, 

OY=d, 

OX=W—d, 
P, the planet, 
8, the satellite. 

8TX= 90° — (i? + ii), 

TXO=im° -~k, 

TXS = 180° ^{p — k), 

T0Z=21Q° — {a—N). 
I am indebted to Professor Frisby for verifying my derivation of the 
formulae (12) and (13), which will be found, I hope, correct, and in a form 
convenient for use. 




